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Abstract 

'Soft' high-energy interactions are clearly important in pp collisions. Indeed, these 
events are dominant by many orders of magnitude, and about 40% are of diffractive ori- 
gin; that is, due to elastic scattering or proton dissociation. Moreover, soft interactions 
unavoidably give an underlying component to the rare 'hard' events, from which we hope 
to extract new physics. Here, we discuss how to quantify this contamination. First we 
present a brief introduction to diffraction. We emphasize the different treatment required 
for proton dissociation into low- and high-mass systems; the former requiring a multi- 
channel eikonal approach, and the latter the computation of triple-Pomeron diagrams 
with multi-Pomeron corrections. Then we give an overview of the Pomeron, and explain 
how the QCD (BFKL-type) Pomeron is the natural object to continue from the 'hard' 
to the 'soft' domain. In this way we can obtain a partonic description of soft interac- 
tions. We introduce the so-called KMR model, based on this partonic approach, which 
includes absorptive multi-Pomeron corrections that become increasingly important as we 
proceed further into the soft domain. This model is able to describe total, elastic and 
proton dissociation data, and to predict the survival probability of large rapidity gaps to 
soft rescattering — in terms of a few physically-motivated parameters. However, more 
differential phenomena, such as single particle p t distributions, can only be satisfactorily 
described if hadronization effects are included. This is achieved by incorporating the KMR 
analytic approach into the SHERPA Monte Carlo framework. It allows a description of 
soft physics and diffraction, together with jet physics, in a coherent, self-consistent way. 
We outline the structure, and show a few results, of this Monte Carlo, which we call 
SHRiMPS, for reasons which will become clear. 
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1 Introduction 



There is no unique definition of diffraction [1]. We may say diffraction is elastic (or quasi- 
elastic) scattering caused, via s- channel unitarity, by the absorption of components of the 
wave functions of the incoming particles. For example at the LHC, there are processes pp 
pp, pX, XX, where neither, one, or both, protons are allowed to dissociate into a system 
X with the quantum numbers of the proton. This definition is only a useful definition for 
quasi-elastic processes, but not for dissociation into high-mass systems. 

An alternative definition is that a diffractive process is characterized by a large rapidity 
gap (LRG), which is caused by t-channel Tomeron' exchange. However this definition is only 
good for very LRG events; otherwise the gap can be due to secondary Reggeon exchange or by 
fluctuations in the hadronization process. 

Let us return to s-channel unitarity: SS^ = I with S = I + iT, so that T - T f = iT^T. For 
fixed impact parameter, 6, if we sandwich the relation between elastic states we obtain 

2ImT e i(>, b) = |T e i(s, b)\ 2 + G inel (s, b) which is satisfied by T e i(s, b) = i(l - e~ n/2 ). (1) 

At high energies the elastic amplitude is dominantly imaginary, so the eikonal/opacity, Q is 
real, and greater than or equal to zero. The total, elastic and inelastic cross sections are easily 
expressed in terms of Q. For example 

~^b~ ~ ^b~^b ~ 2ImT ^-|^l ~ l ~ e ' ( 2 ) 

so e~ n ^ b) is the probability of no inelastic interaction at b. The textbook example of how 
absorption (that is the Gn ne i term in (1)) gives rise to elastic scattering, is scattering on a 
black disc with ImT e i = 1 for b < R] then we have <Ti ne i = 7ri? 2 , which induces a e \ = 7ri? 2 , so 
a tot = 2vri? 2 . 

So much for elastic unitarity. What about proton dissociation? To include dissociation into 
low-mass states we have the Good- Walker formalism. We introduce combinations of £>,£>*, 
(the so-called diffractive eigenstates z, fc, ..) which only undergo 'elastic' scattering. Thus (1) is 
generalized into a multichannel eikonal formalism, sketched in Fig. 1(a). 

What about dissociation into systems of high-mass Ml High M production may be rep- 
resented by Pomeron exchange, giving rise to the triple-Pomeron exchange diagram, shown 
symbolically in Fig. 1(b) with its multi-Pomeron absorptive corrections. 



2 Why study diffraction? 

Well, first, there is its intrinsic interest; the LHC offers the possibility to probe the asymptotic 
behaviour of the pp interaction. Here, we are concerned with a more practical reason. About 
40% of the events at the LHC are due to diffraction — elastic scattering and proton dissociation. 
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Figure 1: High energy pp scattering, showing (a) the multichannel eikonal formalism to allow for 
low-mass proton dissociation, and (b) the triple-Pomeron diagram and multi-Pomeron corrections 
which account for high-mass dissociation. 

Together with other 'soft' interactions, these events dominate by many orders of magnitude. 
Moreover the LHC detectors do not have 47T geometry. Although experimental triggers are able 
to select the rare 'hard' interaction events at the LHC, from which we seek evidence of 'new 
physics', they contain particles from the 'soft' underlying interaction. It is therefore of great 
importance to construct a Monte Carlo, which includes diffraction, and which merges soft and 
hard interactions in a coherent self-consistent way. 

A simple example, illustrating the need of such a Monte Carlo, is related to the measurement 
of the energy scale of jets - additional hadronic activity beyond the standard parton shower 
plus hadronization approach clearly may change measurements of, say, the p± distribution of 
jets. Another example is the potentially informative exclusive production processes at the LHC 
where one can study the quasi-elastic hard subprocess in a very clean experimental environment. 
Unfortunately, in the present LHC experiments we are unable to determine the full kinematics 
of these exclusive events. Such events are selected simply by the existence of a Large Rapidity 
Gap (LRG). Thus we need a reliable MC generator to estimate the mixture of the pure exclusive 
process with the processes where an incoming proton dissociates into either a low-mass or a 
high-mass system. Such a Monte Carlo would also be valuable for interpreting high energy 
cosmic ray data as obtained, for example, in the Auger experiment. 

To construct the sought-after Monte Carlo, which describes soft and hard interactions in a 
unified framework, we first need a partonic model of the Pomeron. 

3 The Pomeron 

Conventionally, 'soft' and 'hard' high-energy pp interactions are described in terms of different 
formalisms. High-energy soft interactions are described by Reggeon Field Theory [2] with 
a phenomenological (soft) Pomeron, whereas for hard interactions we use a QCD partonic 
approach. In perturbative QCD (pQCD), the Pomeron is associated with the BFKL vacuum 
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singularity [3]. However, the two approaches appear to merge naturally into one another. That 
is, the partonic approach seems to extend smoothly into the soft domain. 

The BFKL equation describes the development of the gluon shower as the momentum 
fraction, x ) of the proton carried by the gluon decreases. That is, the evolution parameter is 
ln(l/x), rather than the ln/c^ evolution of the DGLAP equation. Formally, to justify the use 
of pQCD, the BFKL equation should be written for gluons with sufficiently large k t . However, 
it turns out that, after accounting for next-to-leading ln(l/x) corrections and performing an 
all-order resummation of the main higher-order contributions [4], the intercept of the BFKL 
Pomeron depends only weakly on the scale for reasonably small scales. The intercept is found 
to be a P (0) =~ 1.3 over a large interval of smallish k t . Thus the BFKL Pomeron is a natural 
object to continue from the 'hard' domain into the 'soft' region. 



Figure 2: The cascade structure of a gluon ladder. The BFKL or QCD Pomeron is the sum of 
ladder diagrams, each with a different number of rungs. 



In terms of Feynman diagrams, the BFKL or QCD Pomeron may be viewed as a ladder 
diagram built by the exchange of two i-channel (Reggeized) gluons, see the left-hand side of 
Fig. 2. The sequence of parton spitting produces a gluon cascade which develops in ln(l/x) 
space, and which is not strongly ordered in k t) see the right-hand side of Fig. 2. There are 
phenomenological arguments (such as the small slope of the Pomeron trajectory 2 , the success 
of the Additive Quark Model relations, etc.) which indicate that the size of an individual 
Pomeron is relatively small as compared to the size of a proton or pion etc. Thus we may 
regard the cascade as a small-size 'hot-spot' inside the colliding protons [5]. 

At LHC energies the interval of BFKL ln(l/x) evolution is much larger than that for DGLAP 
ln/c 2 evolution. Moreover, the data already give hints that we need contributions not ordered 
in h t , a la BFKL, since typically DGLAP-based Monte Carlos, tuned to describe the Tevatron 
data, overestimate the observed (k t ) and underestimate the mean multiplicity observed at the 
LHC [6, 7]. Further, it is not enough to have only one Pomeron ladder exchanged; we need to 
include multi-Pomeron exchanges. 

Basically, the picture is as follows. In the perturbative domain we have a single bare 'hard' 
Pomeron exchanged with a trajectory «p are c± 1.3 + o^ are £, where c^ are ^0.05 GeV -2 . The 

2 Recall that the slope a' P oc oc Rp om . 
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transition to the soft region is accompanied by increasing absorptive multi-Pomeron effects as 
we go to smaller k t , such that a soft amplitude may be approximated by an effective linear 
trajectory afp ~ 1.08 + 0.25i in the limited energy range up to Tevatron energies [8]. This 
smooth transition of the Pomeron from hard to soft, that is from bare to effective, is well 
illustrated, for example, by the behaviour of the data for vector meson (V = p, cj, 0, J/i/j) 
production at HERA, 7*(Q 2 ) +p -> V(M) + p. As Q 2 + M 2 decreases from about 50 GeV 2 
towards zero, the s and t dependence of the data reveal the trend that the effective intercept, 
afF(O), decreases from about 1.3 to 1.1, and the slope, cV P , increases from about zero to about 
0.2 GeV- 2 . 



4 Partonic structure of the Pomeron — the KMR model 

How can we implement the partonic model of the Pomeron in practice? This is achieved [9] 3 
from the partonic ladder structure of the Pomeron, Q(y, k t) 6), generated by BFKL-like evolution 
in rapidity. y ) with a simplified form of the kernel, but which incorporates the main features 
of BFKL: diffusion in ln/c 2 and A = ap are (0) — 1 ~ 0.3. At each step of the evolution in y, the 
value of the partonic k t and the impact parameter, 6, can be changed. 

It is important to note that, in comparison with Reggeon Field Theory, Q now depends on 

an extra variable, k t) that is the transverse momentum of the intermediate parton, in addition 

— * 

to the usual two 'soft' Regge variables: y = ln(l/x) and the impact parameter, 6, which is 
conjugated to the transverse momentum, Q t , transferred through the entire ladder. 

After including rescattering of the intermediate partons, we actually have two coupled evo- 
lution equations. One evolving up from the target k at y — 0, and one evolving down 4 from 
the beam i at y' — Yj — y = with Yj ln(s/k' 2 ), 

exp(-X[Q k (y) + W)]/2) K{k u k[) Q k (y). (3) 





f d 2 k't 


dy 






f d 2 K 


dy' 


1 *k? 



exp(-A[^( ? / / ) + Q k (y)]/2) K(k u k' t ) ^(j/), (4) 

where, for clarity, we have suppressed 5 the k t labels of the JTs. The inclusion of the k t 
dependence is crucial for the transition from the hard to the soft domain. The absorptive 



3 A two-channel eikonal is used, z, k = 1, 2. 

Yj accounts for the fact that au ±ci± 6 c± ro t 



4 The value of Yj accounts for the fact that at larger k[ a smaller rapidity interval is available for the evolution. 
5 In general, the kernel K depends on the difference b — b' . However, the b dependence is proportional to 
the slope a' of the bare Pomeron trajectory. Asymptotically the BFKL approach predicts a' — » 0; and indeed, 
analyses of soft data find it to be very small. Thus, for simplicity, we have not shown this variation in eqs. (3, 
4). Then the only b dependence of Q comes from the starting distribution of the evolution, and not from the b 
dependence of K. 
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(exponential) factors 6 in the equations embody the result that there is less screening for larger 
k t . The coupled evolution equations may be solved iteratively to give Q(y,k t ,b), for a whole 
range of fixed values of b. 

Recall that the absorptive factors, exp(— in (3, 4) lead to a strong suppression 
of the low k t domain, which introduces an effective, dynamically generated, infrared cutoff 
kt > fcmin, whose value increases with energy. Therefore the final high energy results are not 
too sensitive to the value of the original cutoff that is included artificially in the propagators 
(l/k 2 — >• l/(k 2 + Qq)) in order to protect the infrared singularity. 

In principle, knowledge of Qk(y,k t ,b) allows a good description of all soft and semi-hard 
high-energy pp data, such as cr to t, da e \/dt, dasB / dtdM 2 , ..; and reliable estimates of the survival 
factors of large rapidity gaps [11]; and even semi-quantitative estimates of PDFs and diffractive 
PDFs at low x and low scales. Moreover, such a model [9] has only a few physically moti- 
vated parameters. Specifically, the parameters are the bare Pomeron intercept A = op(0) — 1 
(expected to be about 0.35) and slope a' P (expected to be small); d which controls the diffu- 
sion in Ink 2 ] A; the initial gluon density, TV; and the parameters which specify the diffractive 
eigenstates. 

Examples of topical predictions of the model are that the rapidity plateau increases with 
collider energy as da/dy ~ s 0-2 , just like the LHC data from 0.9 to 7 TeV; or that the gap 
survival for double diffractive (exclusive) SM Higgs production, pp p + h + p at 14 TeV the 
LHC 7 is about 0.02. Exclusive (low multiplicity) processes, like da e \/dt and pp —> p + h + p ) 
may be calculated analytically, whereas in inclusive cases we need a Monte Carlo to account for 
hadronization effects. In particular, we need to allow for hadronization to describe the ATLAS 
measurements, Fig. 3, of da/d(Arj) for Ar] < 5; although for larger gaps the triple- Pomeron 
diagram with absorption corrections will suffice, see the discussion in [16]. 

5 The SHRiMPS Monte Carlo 

Most of the existing general purpose Monte Carlo event generators (Pythia, Herwig, Her- 
wig++, Sherpa), [17], use collinear factorisation and standard 2—^2 QCD matrix elements 

6 Since we are dealing with the amplitude fJ, and not with the cross section, we use here exp(— AO/2) and not 
exp(— XQ). The parameter A determines the strength of the triple- (and multi-) Pomeron couplings in terms of 
the Pomeron-proton coupling. In practice, a slightly different form of the absorptive factor [9] is preferred by 
the data, which is consistent with the AGK cutting rules [10]. 

7 If the outgoing protons in this process are measured far from the interaction point then the mass of the 
SM Higgs particle can be determined, via missing mass, to an accuracy of about 1 GeV. If is less than 
about 140 GeV, this process offers a chance to study the h — » bb decay, since the QCD bb background is heavily 
suppressed by a J z = selection rule [12]. Such exclusive events are very clean, but a ~ 2 fb [13], and events are 
contaminated by 'pile- up'. With moderate pile- up, very accurate event timing can select the exclusive events. 
An advantage of such a process is that, in the case of MSSM, the h and H Higgs of may be both observable, 
and the cross sections may be enhanced. 
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Figure 3: The inelastic cross section differential in rapidity gap size Arj for particles with p T > 200 
MeV. The data are from [14]. The triple-Pomeron diagram, with multi-Pomeron corrections, can 
describe the high Ar] data [15], but hadronization effects are evident at smaller Ar]. 



supplemented with a factorised spatial dependence as the basis for their eikonals, Q. In order to 
avoid potential problems in the infrared region, however, these matrix elements are either cut 
at some low value of minimal k t of the outgoing partons, or the soft behaviour of the i-channel 
and as are modified (\t\ —> \t\ + q^) to render the result finite. In the former case, when the 
'hard' eikonal is cut at fc , it is further supplemented with a soft one, Q = £\ard + fi so ftj to 
account for the contribution from the (infrared) region k t < k^. In any case, the simulation 
typically then uses total cross sections etc. as input parameters and concentrates on producing 
the correct particle distributions. 

Here, we seek a Monte Carlo that describes all aspects of minimium bias events — total, 
differential elastic cross sections, proton dissociation, diffraction, jet production etc. — in a 
unified framework, capable of modelling exclusive final states. In order to obtain such a Monte 
Carlo, which includes a description of soft physics and diffraction as well as jet physics in a 
self-consistent way, the KMR model has been incorporated into the SHERPA framework [18]. 
Recall that the KMR model is based on the bare QCD Pomeron, with absorptive multi-Pomeron 
rescattering corrections, that increase as we continue from the hard to the soft domain. For 
this reason we call the Monte Carlo 'SHRiMPS', standing for Soft Hard Reactions involving 
Multi Pomeron Scattering. 

Let us outline the construction of the Monte Carlo. First, we solve the coupled evolution 
equations in rapidity y to generate Qik(y, k t , 6), having specified the boundary conditions of the 
diffractive (Good- Walker) eigenstates. These eigenstates specify the elastic and quasi-elastic 
scattering. For a detailed simulation of the inelastic state we must select the number of ladders 
to be exchanged. We take the number to be distributed according to a Poisson distribution, 
with parameter 0^(6). Having decided on the number of (primary) ladders to be exchanged, 
the incoming protons are dissolved into a valence quark, a (non-interacting) valence diquark and 
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gluons. A random pair of partons from these beam dissociations is chosen to exchange the next 
ladder. Gluon emmisions are generated in the ladder according to a Markov chain, ordered 
in rapidity, with a pseudo-Sudakov form factor. The i-channel propagators in the ladder are 
reggeized gluons that come in different colour states; only singlets and octets are considered. 
For each t-channel propagator we select the probability of a colour singlet, 

Vl = (1 _ e -«V2 )2) (5) 

where SQ is known from the solution of (3, 4). We take V$ = 1 — V\. Each gluon emission 
leads to two new propagators: allowed combinations are V{P%^V^P\ and V%P%. The formation 
of singlet propagators gives rise to rapidity gaps being associated with 'elastic rescattering'. 
Simultaneously the rescattering may lead to the inelastic interaction of secondaries, produc- 
ing new ladders with Poisson probabilities e~ sn (5Q) n jn\. The single gluon emission process is 
iterated until the active interval becomes a colour singlet or no further emissions are kinemat- 
ically allowed in the rapidity interval. Finally we implement the 'usual' parton shower, plus 
hadronization, plus hadron decays, plus QED to produce the final scatter of observed particles. 



Toward £ p ± density vs. /> trk \ y/s = 7 TeV Away £ p ± density vs. p* kl , y/s = 7 TeV 




Transverse T y 1 density vs. y/s = 7 TeV 




Figure 4: The values of J2pt in the toward, away and transverse regions to the direction of a 
leading jet of given p T (jet). The data are from [19]. 

The Monte Carlo, which has relatively few parameters (several of which are physically 
motivated), is tuned to an extensive range of 'soft' high energy data. Two examples of the 
preliminary description of data are shown in Figs. 3 and 4. 
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